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Abstract 

The nonrelativistic Schrodinger equation for motion of a structureless 
particle in four-dimensional space-time entails a well-known expression 
for the conserved four- vector field of local probability density and current 
that are associated with a quantum state solution to the equation. Un- 
der the physical assumption that each spatial, as well as the temporal, 
component of this current is observable, the position in time becomes an 
operator and an observable in that the weighted average value of the time 
of the particle's crossing of a complete hyperplane can be simply defined: 
the theory predicts, and experiment is presumed to be able to observe, 
the integral over the hyperplane of the normal component of probability 
current, weighted by the time coordinate. In conventional formulations 
the hyperplane is always spacelike, i.e., is a time=constant hyperplane in 
Galilean relativity, and the result is then trivial. A nontrivial result is 
obtained if the plane is not of this type. When the space-time coordinates 
are (t,x,y,z), the paper analyzes in detail the case that the hyperplane 
is of the type z=constant. Particles can cross such a hyperplane in ei- 
ther direction, so it proves convenient to introduce an indefinite metric, 
and correspondingly a sesquilinear inner product with non-Hilbert space 
structure, for the space of quantum states on such a surface. Since the 
metric is indefinite, an uncertainty principle involving the dispersion of 
the crossing time and the dispersion of its conjugate momentum does not 
appear to be derivable from the theory. A detailed formalism for com- 
puting average crossing times on a z=constant hyperplane, and average 
dwell times and delay times for a zone of interaction contained between a 
pair of z=constant hyperplanes, is presented. 
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1 Introduction 



Within a few years after the discovery of quantum mechanics a consensus formed 
(von Neumann [I], p. 188 and p. 354, Pauli 0, P- 140, footnote, and 0], p. 
63, footnote) to the effect that, in contrast to spatial positions, and therefore in 
conflict with special relativity, the temporal position t is necessarily a c-number, 
or parameter, with no generic operator status being mathematically feasible. In 
the decades intervening since the publication of the original versions of the two 
cited treatises in 1932 and 1933, respectively, the prohibition on specifying the 
time as a dynamical variable has been widely upheld as part of the standard 
doctrine of quantum mechanics — see, e. g., Peres [Sj Chap. 12-7, Onmes 6 p. 57, 
and Sakurai 7; p. 68. In recent decades interest in this subject has intensified, 
due in part to applications of tunneling phemomena in semiconductors ([5], 
P] Ch. 3.2.3, ^Hli |H| - |12| . |13p. and a substantial set of results has been 
published that introduce formalisms that argue for, or against, various quantum- 
mechanical definitions of time, including arrival times, tunneling times, dwell 
(alias sojourn) times, and delay times. For definitions, reviews, and citations, 
see Refs. Q3], US], P2|, ESI, d, and a collection of articles in |20) . 

This paper claims to advance a general operator for time as a quantum- 
mechanical observable in the context of the nonrelativistic Schrodinger equa- 
tion. A natural preliminary question is, how does the construction of such an 
observable time square with the putative result that no such operator exists 
in general? I believe that the nonexistence proof rests on an improper substi- 
tution of dynamics for kinematics, and that there is a straightforward inter- 
play between classical and quantum mechanics on this question. We infer from 
Pauli 3 in a footnote on p. 140, that the search for a classical observable for 
time reduces to the following problem: Let {q 1 ,Pi, ■ ■ ■ , q N ,Pn} be the positions 
and momenta in a generic classical Hamiltonian dynamics, with Hamiltonian 
-ff(i, g 1 ,pi, . . . ,q N ,pn)- We define the ordinary Poisson bracket {/i,/2} P b of 
two functions f\,2{t, q 1 ,Pi, ■ ■ ■ , q N ,Pn) to be 

U f\ ^r( d hdf-2 dhdf 2 \ 

{/i,/ 2 } Pb - QwwrWiWr { ) 

Then we want to find a function T(t, q 1 ,p±, . . . , q N ,pn) such that 

{T,ff} pb = 1. (2) 

Such a function T would be the classical limit of a quantum-mechanical Her- 
mitcan operator T qm that is invoked by Pauli, but demonstrated by him not to 
exist in general, such that the commutator satisfies 

Pqm, ^qoi] = ^lqm- (3) 

Pauli's theorem has prevailed over the intervening decades, and attempts to de- 
fine an observable time in connection with solutions to the Schrodinger equation 
have had recourse to many alternate approaches — see the reviews cited above. 
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It is easy to see in a related context that Pauli's result, although correct in 
its limited context, addresses the wrong question. Let us define the augmented 
Poisson bracket {Fx, i^japb of two functions Fi j2 (£,Pt, q l ,p%,..., q N ,pn) to be 

dFi dF 2 dFi dF 2 
{Fi,F 2 } apb = ——-—— + {Fl ,F 2 } ph . (4) 

The function T> c \, defined as 

V cl = Pt + H, (5) 

generates the time dependence of quantities as F, in that 

dF 

= {^cl}apb, (6) 

subject to the dynamical constraint that on classical paths we require 

£>d = (7) 

(see |2I], Ch. VI.10). We now infer that i£F = t, then dF/dt = {i,X> c i}a P b = 1, 
whatever be the Hamiltonian H . In quantum mechanics, with suitable attention 
to operator ordering, we substitute 

h d i h d 

V qm = m + H(W, 1W ,.:), (8) 

and obtain the time-dependent Schrodinger equation for ip(t,q l , . . . ,q N ) as a 
dynamical constraint on ip, that is, 

Pqm^ = 0. (9) 

If ip is a solution to the Schrodinger equation we can substitute Hip for ihdip/dt, 
but not otherwise. In particular, if tp is a nontrivial solution to 10, tip is not a 
solution, so that 

d 

ifi—ttp ^ Htip. (10) 

Therefore, trying to satisfy the operator commutation rule J2J is not relevant 
to the problem of finding an operator for the time, even when the operand is 
a solution to the Schrodinger equation. However, the general definition for the 
time-derivative of quantum-mechanical operators 

dF, 



_jp = (01)-^^} (ii) 

does make sense no matter what the state function, dynamically constrained or 
not, operated on by the rhs: in fact, if F qrn = t, whatever be the Hamiltonian, 
its time derivative operator is the unit operator. 

Let the classical Hamiltonian for a particle in four-dimensional space-time 

be 

H = ^-{p 2 x +p 2 y +p 2 z )+V(t,x,y,z). (12) 
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Then (JJJ can be solved for p z to yield 

^dit, Pt,x,p x ,y,p y ,z,p z ) = p z T [-2m(p t + V) ~pl~pl} 1/2 = 0. (13) 

The quantity T>' ±cl , considered as a function of eight variables prior to its role as 
a constraint, generates the equations of motion for the evolution of the system 
in the indirection. The same classical dynamics is obtained if we return to 
the original Lagrangian formalism, take t as a dependent variable and z as the 
independent variable, and procede to the Hamiltonian formalism. We remark 
that the quantum operator hH zcv of l|24|l can be construed as a matrix square 
root of the quantum operator that arises from a term in l|13[l . 

In this paper we shall consider only the nonrelativistic form of quantum 
mechanics, restricted to the problem of determining the wave function of a 
massive, structureless particle in Galilean four-dimensional space-time in the 
presence of a given space- and time-dependent potential energy distribution. 
The physical hypothesis that underlies the theory herein is that not just the 
time component, but also the spatial components, of the conserved four-current 
density of probability flow are observables. This hypothesis does not seem to 
have been made explicitly, or its consequences studied, heretofore. A sketch of 
the mathematics introduced to elaborate this physical assumption follows. 

We shall treat quantum mechanics as a boundary-value problem for the 
Schrodinger equation. The specified boundary values on a simple domain in 
space-time will be regarded as input to the problem, and the derived interior 
values and complementary boundary values will be considered as the overall 
output. An intermediate objective will be to formulate a theory of spatial 
evolution of wave functions. As is usual in applications, we shall choose one 
spatial coordinate as the evolution coordinate, and expand the wave function in 
terms of conveniently simple orthogonal functions of the transverse coordinates, 
which include the time. The time and the two transverse spatial coordinates 
therefore will appear naturally as operator/observables in the space of such 
functions, analogous to the role of the three spatial positions when the time is 
taken as the evolution coordinate. 

Suppose that we want to solve the Schrodinger equation in a space-time box 
surrounded by two i=constant walls and a spatial boundary. Since the equa- 
tion is of first order in time and of second order in the spatial coordinates, it 
is, mathematically speaking, necessary and sufficient to supply wave function 
values on the earlier i=constant surface, and a suitable combination of wave 
function and normal-derivative values on the spatial boundary surface, to infer 
that an interior solution exists and is unique, as discussed in [22], Ch. 5, §3. 
(Our mathematics differs from Friedman's in that his equation is the heat equa- 
tion rather than the Schrodinger equation, and we shall administer nonlocal 
boundary conditions, which distinguish input from output signals, on the spa- 
tial boundary.) Conventional time-dependent quantum mechanics for the most 
part deals with specifying initial, or (but not and) occasionally final, values on 
a t=constant surface and simple (often, zero) values on the spatial boundary, 
which can be partly or wholly at infinity. Nontrivial spatial boundary values, 
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as incoming wave amplitudes in a scattering problem, are conventionally spec- 
ified only in the context of the time-independent Schrodinger equation. In the 
present work we shall generalize the latter problem by considering general time- 
dependent, as well as space-dependent, input values on the spatial boundary, in 
the presence of explicit time dependence in the potential energy function in the 
differential equation itself. We shall hereinafter denote these cases of boundary 
value problems as Type I and Type II, respectively. These correspond, roughly 
and respectively, to the first and second initial-boundary value problems ana- 
lyzed in [22], Chs. 3 and 5. 

In the first problem, the wave function evolves in time from given initial 
values, with time-independent spatial boundary values. In the second problem, 
we shall consider that the wave function evolves with respect to the z-coordinatc, 
such that the interior domain corresponds to a finite interval in the chosen 
coordinate z. Since the differential equation is of second order, determining the 
evolution of a wave function in a spatial direction is generally a more difficult 
task of analysis in both the mathematical and physical senses than one for 
its evolution in time. We summarize the derivation to be carried out below 
in terms of the following ten observations, steps, or results: (i) the space of 
states on any given z = z\ hyperplane has a natural doubled structure in that 
it comprises the direct sum of the values and of the ^-derivatives of the usual 
space of wave functions ij)(t, x,y, z) at z = Z\\ (ii) the Hamiltonian is a 2 x 2 
matrix of operators that is derived from the ordinary Schrodinger equation; (iii) 
the familiar expression for the probability current density in the z-direction is 
used to infer the definition of a metric operator in the space of states, where 
now inner products include an integral over t as well as over x and y; (iv) the 
metric so derived is indefinite, and the Hamiltonian is self-adjoint with respect 
to the metric (synonymously, pseudo-Hermitean); (v) the norm being indefinite, 
we shall sometimes use the term "particle presence" to denote the unit operator, 
the expectation value of which is the above-mentioned norm; (vi) apart from 
modifications needed for closed channels, the formalism can be established so 
that waves traveling in the +z direction have positive norm, and waves traveling 
in the —z direction have negative norm with respect to the metric; (vii) while for 
open-channel modes the direction of travel and of propagation will coincide in 
the large, it is convenient to define these categories differently for closed-channel 
modes; (viii) the input and output at either end of a finite spatial interval [zi, z^\ 
are taken to comprise, respectively, the superposition of waves propagating into, 
and the superposition of waves propagating out of, the interval at the initial 
point z\ and at the final point z-i (this means that there will be only outward 
propagating scattered waves from a zone of interaction); (ix) orthonormal sets 
of input or of output states, transition amplitudes, and probabilities are then 
computed using what amounts to a Hilbert space inner product, which is derived 
from the indefinite metric, but depends on the wave function and its z-derivative 
at both Z\ and z^\ (x) the dynamics yields a mapping of open-channel input into 
open-channel output that is unitary. 

The fact that a pseudo-Hermitean Hamiltonian describes the spatial evolu- 
tion of a physical system's wave function has another concomitant: the Hamil- 
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toman can have, as well as real eigenvalues, nonrcal eigenvalues that occur in 
complex conjugate pairs (Gohberg, et al., (231 P- 23, Proposition 2.4). Each 
such pair is associated with the two wave function solutions (one rising, the 
other falling exponentially) for a closed channel, or classically inacessible re- 
gion for the system when it is in an associated quantum state. We shall argue 
that it is natural to define the direction of propagation (but not of travel of 
the particle that the wave represents — see the discussion in Sec. 4) of such a 
wave as the direction in which it decreases exponentially. The simple exponen- 
tial states in such a pair each have zero norm and, with proper normalization, 
unit overlap, which complicates the formalism. A further complication results 
from the circumstance that a degenerate eigenvalue of a Hamiltonian requires 
special treatment when the Hamiltonian cannot be diagonalizcd by a similarity 
transformation, leading to the appearance of so-called "n-pole ghost" quantum 

states, for n = 2, 3, There is more discussion on these problems below. 

Formalisms for the spatial evolution of a wave function were proposed by 
Kijowski |21] and by Piron [55], and their work was discussed by Mielnik |2(i|. 
These two approaches differ substantially from each other and from the for- 
malism introduced herein, as will be discussed following Eq. I|24|) and in Sec. 
4. 

The quantum mechanics describing evolution of a wave function in both 
directions across a spatial interval is to an extent patterned after the author's 
previous work |27| on a quantum dynamics that encompasses joint bidirectional 
evolution of a quantum state between two temporal walls. 

The remaining sections are organized as follows: In Sec. 2, we shall formu- 
late expressions for the four-current density associated with a physical quantity, 
and for the local space-time density for creating or destroying that quantity in 
a quantum-mechanical system. We shall also show how to prescribe physically 
motivated boundary conditions so that the Schrodinger equation can be solved 
in a semi-infinite (finite in the z-direction, infinite in the t, x, y-directions) box. 
In Section 3 we shall propose a formalism for computing the average temporal 
position (i.e., crossing time) of the particle at both spatial walls of the box, 
given the spatial input and given the 5*-matrix deriving from a general interac- 
tion potential energy in the box's interior. These results will then be used to 
compute formulas for dwell and delay times for the particle remaining within, 
reflecting from, or transmitted across, the box. Section 4 concludes the paper 
with a discussion of the present formalism and of previous work on the subject. 
The matters of an uncertainty principle involving the time and its conjugate 
momentum, and of an extended theory of measurement, are merely touched 
on there, as attempts by the author to develop these constructs have not been 
successful. 

2 Quantum-mechanical formalism 

In this section we shall set up the theory that forms the "floor" of the present 
work. Rather than attempt to make the formalism highly general, we shall 
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develop the argument in a particular context: the wave function solution of 
Schrodinger's equation for a particle moving in the interior of a certain simple 
box of four-dimensional space-time. In particular, we propose a formalism and 
an interpretation that incorporate the wave function into an expression for the 
space-time four-vector "flow" density of a physical quantity, which quantity 
corresponds to a certain linear operator in the function space of fully time- and 
space-dependent wave functions. We shall argue that is is natural to regard 
the four-divergence of the flow as the local density of creation and destruction 
of that quantity at a point in space-time for the physical system in that time- 
dependent quantum state. Either the volume integral of the divergence, or the 
surface integral of the normal component of the flow vector density, therefore 
represents the total amount of that quantity generated inside the space-time 
box. If that quantity is the time t, this integral plausibly represents the average 
so-called dwell time of the particle in the given box, given that the wave function 
is properly normalized. 

Let B\ and B2 be the following open boxes in space-time: 

B\ = {(t, x, y, z)\t\ < t < ti, —00 < x < 00, —00 < y < 00, —00 < z < 00}. 

(14a) 

$2 = {(£) x, y, z)\ — 00 < t < 00, —00 < x < 00, —00 < y < 00, z\ < z < Z2}. 

(14b) 

The Schrodinger equation for ip(t, x, y, z) for a particle of mass m can be derived 
from a variational principle for an action A, as given in Schiff ([22], p. 499), but 
modified to make it real: 



A = dt dx dy dz 



B 



ih difj ihdifj* h 2 

at -"2^-2^^ v &*>v>*)* 



(15) 

The equations of motion are to be obtained by keeping the boundary values of 
tp and tp* fixed, and pretending that in the interior region ip and -0* can be var- 
ied independently and "arbitrarily". The action is stationary when ip(t, x, y, z) 
satisfies 

lt-^ + V ^^ = °> (16) 
and ip(t,x,y,z)* satisfies the complex conjugate equation, for all (t,x,y,z) € 
Bx or B 2 - 

Henceforth when we say "solution" , we shall mean a function defined over 
the entire box such that it satisfies equation (|16|l everywhere in B\ or £>2- The 
linear operators representing physical quantities will normally carry a solution 
into another space- and time-dependent function that is not a solution, so in 
effect we shall deal with the more general vector space of well-behaved, complex- 
valued functions of space and time that need not be solutions of the Schrodinger 
equation. 

A standard problem in conventional quantum mechanics arises if we con- 
strain tp(t, x, y, z) to be zero on the infinite parts of the spatial boundary of the 
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box B\, that is 



and require that 



tp(t,x,y,z) -> 0, if |x| + \y\ + \z\ 



tp(t,x,y, z) — > u(x,y, z), as t 



(17) 



(18) 



where u(x,y,z) is some given complex- valued function on the earlier temporal 
boundary of the box. As is well known, the interior values of ip(t,x,y,z), and 
the limiting values on the temporally later boundary of the box at t = t2, are 
all uniquely determined by the differential equation and these input boundary 
conditions. The derived values are all output in a sense, but we shall often mean 
by output just the subset of boundary values that were not given as input, in the 
present case ipfe, x, y, z) with — oo < x < oo, — oo < y < oo, and -co < z < oo. 

Now let us consider a problem in box £>2 such that certain information about 
the limiting values of ip(t, x, y, z) and dip/dz(t, x, y, z) on the two walls z = Z\ 
and z — Z2 is given as input, while the wave function is supposed to tend to 
zero as t and/or x and/or y tend to ±oo. We want to specify just enough input 
information so that a solution satisfying the input boundary conditions exists 
and is unique. In order to accomplish this, we need to do some preliminary 
work. We shall not keep to a mathematically rigorous derivation, but appeal to 
plausibility arguments at most steps. 

We now convert the above variational principle to Hamiltonian form using 
the methods of Goldstein (UHli Chap. 12-4), with the proviso that it is the 
spatial parameter z, rather than t, that is taken as the evolution coordinate for 
the wave function. The quantity in square brackets in l|15fl is the Lagrangian 
density C The canonical field momenta are 



dC 



Pi>* = 



dip' 



h 2 dip* 
2m dz 

h 2 dj)_ 

2m dz 



The action functional becomes 



A 



dt dx dy dz 



B-2 



dtp dip* 

P^-g- + -g~Pi>' -rt{lp,Pi,,1p 



(19a) 
(19b) 

(20) 



where the Hamiltonian density is 
2m 



H 



n 2 



ih dip ihdip* , h 2 



dip* dtp dip* dtp 
dx dx dy dy 



-iP*Vip. (21) 



The equations of motion obtained by varying ip* and are a coupled set of 
linear equations; a complex conjugate set is obtained by varying ip and p^* . We 
write the former equations in 2 x 2 matrix-operator form as follows: We first 
define 

ip(t,x,y,z) 



(22) 
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then the equations of motion can be written 



19* 

-■5- = H zev ^, 
i oz 



(23) 



where the Hamiltonian H zev (the subscript "zev" stands for "z-evolution" ) is 



in 3 + *L(gs 

ot 2m V ot^ 



o 



V(t, x, j/, 



2/n i 

h 2 




(24) 



We note that Piron obtained a Schrodinger equation for a wave func- 
tion's evolution along the spatial coordinate x, but Piron's wave function has 
one component, and the Hamiltonian is the operator derived from the classical 
quantity that generates dynamical motion along a the x-axis. Piron thereupon 
obtained a general expression for the evolution in x of the average temporal 
position of a particle in one space dimension, but did not develop the theory 
further. 

In the ordinary quantum mechanics derivable from the variational principle 
Eq. 115(1 , the z-component of the conserved probability four-current density is 
(ESI, P- 27) 

2im V dz dz 

In the present language this expression takes the form 



J 3 (t,x,y,z) 



(25) 



where M is the 2x2 matrix 



M = 



i 
-i 



(26) 



(27) 



Note that the matrix M is Hermitean, has unit square, and has eigenvalues ±1, 
so that it can engender an indefinite metric. By inference, we make a guess for 
an inner product law for two z-propagating states: 



(*i(2);* a (*)) = h- 



dtdxdy [*i(^x,y,z) t M* 2 (i,a;,2;,z)]. (28) 



Note that this formula has the appropriate physical dimensions, in that if 
have the usual dimension length~ 3 ^ 2 , then the above inner product is dimen- 
sionless. 

We shall now argue that the above ingredients can be made into a theory 
of spatial evolution of a Schrodinger wave function. We shall work with the 
case of z-evolution of a wave function in four-dimensional space-time, but gen- 
eralizations to other cases, as radial or reaction coordinates (see jHO]) for the 
(3iV + l)-dimensional space-time involved in an TV-particle wave function, are 
formally straighforward. 
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Let S be the space of functions of type Eq. with some appropriate 

boundary conditions. Wc define the M-adjoint of a linear operator W acting 
on this space as that unique operator W* such that 

(W** x ;*2) = (*i;W* 2 ) (29) 

for all $1,^2 £ 5; in 2 x 2 matrix form, with as the ordinary Hermitean 
conjugate, we have 

W t = MW^M. (30) 

If 

W i = W (31) 
we call W pseudo-Hermitean, and if 

Wi = W- 1 (32) 

we call W pseudo-unitary. 

If W is pseudo-Hermitean and the state \& is suitably normalized, we want 
to make the plausible specification that the (necessarily real) number ("J; W^) 
is the expectation value of W in the state We argue in favor of this axiom 
as follows: Let x° — t, x 1 = x, and so on. Suppose that in the conventional 
Schrodinger formalism, lu is some physical quantity, such as the time t, the spa- 
tial positions x, y, z, or the "particle presence" 1 (we denote operators standing 
for physical parameters with a "hacek" accent over the symbols). We take as a 
physical axiom that the four-vector "flow" density jj^'(t, x, y, z), fj, = 0, 1, 2, 3, 
of lu is 

Jo^CMjj/jz) = ip(t,x,y,z)*ujip(t,x,y,z) (33a) 

_ ^t* 
dx 



-^Zk (*' x ' ^ z ) ) > Ior fc = 1) 2 , 3 - ( 33b ) 



The above expressions need symmetrization if u) contains derivative operators, 
e.g., lu = X-! m = xi + t(iH/m)d/dx : ' , one of the components of the initial center- 
of-mass position for a free particle. We compute the four-divergence of the above 
vector field, assuming that ip is a solution to Eq. iJTSj l: 



h / 9w dip dtp* du 
2im \ dx k dx k dx k dx k 



(34) 



This divergence can be construed to be the local density of creation or de- 
struction of the quantity lu by the system in the state ip(t, x, y, z). If the diver- 
gence is zero, as for the case lu = 1, the associated quantity is not being created 
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or destroyed and is both globally and locally conserved. If u> = t, we find that 



3 dJ (i) 



^=0 



(35) 



Hence, the so-called "dwell" time td of the particle in a box B with the given 
input is the space-time integral of the density of creation of time over the box, 
that is, 



td 



dtdxdydz\ip(t, x, y, z)\ 



(36) 



The latter result reproduces a formula given in Ref. ^2] Eq. (2.2), Ref. Eq. 
(2.67), and Ref. 32 Eq. (14). By the divergence theorem we can convert the 
volume integral to a surface integral, so that we have either 



dx dy dz Jq (t, x, y, z) 
dt dx dy Jg*' (t, x, y, z) 



t=ti 



(37a) 
(37b) 



Hence if we have a boundary problem of Type I, such that ip(t, x, y, z) is zero on 
the spatial walls, and ip has the usual conserved unit norm on the i=constant 
walls, we find, with Eq. 133a(l 



' D 



(38) 



If we have a boundary value problem of Type II — we shall discuss later how to 
normalize ip in that case — so that tp(t, x,y, z) tends to zero as \t\ or |x| or \y\ 
becomes large, and using Eq. I|33bl) then 



il 



I I I dtdxdy[ip(t,x,y,z)*t^(t,x,y 



lira 



dib* , . , . . 

-g-{t,x,y,z)tip{t,x,y,z) 



dz 



(39) 



The alternate forms of the dwell time given in Eqs. (|36|l and l|39ll were previ- 
ously obtained by Jaworski and Wardlaw and applied in a series of papers f |33| 
Eqs. (4.2) and (Al), PU, G2|). Given that we compute x, y, z) from 
ip(t,x,y,z) by Eq. I|22|). and that I2 is the 2x2 unit matrix, then the opera- 
tor i/ 2 is the time operator in the z-evolution formalism, and we have, in the 
notation of Eq. 



-h 1 = ^{z)-{tI 2 )^{z))\ z - = Z. 



(40) 



The above results suggest that for a pseudo-Hermitean operator W in the 
space of ^-solutions, and for a boundary value problem of Type II, we should 
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define the expectation value (W)y( z ) of W in the state ^>(t, x, y, z) at a chosen 



z as 

(W}* (z) = (*(z);W*(z)), (41) 

as was proposed earlier in this section. The value (W)$( z ) therefore (in Type 
II problems) specifies the average net flow of W across the given z=constant 
surface. The difference of the expectation values of W computed at z = Z2 
and z = z\ is therefore the net flow of W out of the box, in other words is, on 
average, the total amount of W "created" by the system in the box. 

We want now to define input and output on the spatial walls of the box. We 
define a complete, orthonormal basis <t>(k t ,k x ,k y )(ti x -,y) f° r au {t>x,y) G K 3 as 
follows: 

<t>(k t ,k x ,k y ){t>x,y) = (2tt)~ 3/2 exp(~ik t t + ik x x + ik y y), (42) 

where kt, k x , and k y each range independently from — oo to +00. Although the 
physical dimension of k t differs from that of k x and k v , it is convenient to use 
three-vector notation k = (kt,k x ,k y ) and call the three-volume element d 3 k = 
dk t dk x dk y . The negative sign before k t in the exponent in Eq. l|4^|) is chosen so 
that positive k t corresponds to positive energy; the conjugate momentum to t 
is p t <-> (h/i)d/dt <-> —hk t . 

In an expansion of a wave function ^(t, x, y, z) in the above basis functions, 
we will encounter certain quantities repeatedly, so we now define simplified 
notation for them: Let £ take cither value F or B, which stand for forward and 
backward propagation along z, respectively. We also take 

-(C) ^ | + 1 iK = F ' (43) 

If (2mkt/h) > (k x + ky) (called an open channel), we define 

k z (k) = [2mk t /h-k 2 x -k 2 y ] 1/2 , (44) 
and if (2mkt/ti) < (k x + ky) (called a closed channel), we define 

K,(k) = [-2mk t /h + k 2 x + kl] 1/2 . (45) 

We shall normally just use k z and k z without explicitly citing their arguments, 
except that primed, double primed, and triple primed arguments will be de- 
noted, respectively, by k' z , k", and k"' , and similarly for k z . 
Let a wave function have the expansion in basis functions 

p 

*{t,x,y,z) = III d 3 kJ2ft(k)Mt,x,y)X C (k;z), (46) 

J J JM 3 ]~~ B 

where the /^(k) are the expansion amplitudes, and where the A^(k; z) are nor- 
malized solutions for forward or backward motion along z, which we construct 
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as follows: Substituting Eq. (|4(j|l into Eq. I|23|l . we find that 



1 dX^ 



i dz 

where, for V a function of z alone, 
-ffzcv(k; z) 



(k;z) = H zev (k;z)Xt(k;z), 







2mi/h 2 
2mV(z)/h 2 } 



(h 2 /2mi)[2mk t /h- k 2 x - k 2 y 
When V(z) = 0, and for open channels, we obtain the solutions 



[™/(te 2 )] 1 /2 e xp[a(C)^^] 
-a(C) (i/2) (S 3 k z /m) V 2 exp [a(C)*fc, z] 



the corresponding inner products are independent of z: 

/T 1 X c '(k;z) t JlO' c (k;z) = <5 c ' c cr(C)- 



(47) 



(48) 



(49) 



(50) 



Note, however, that these solutions do not satisfy the Cauchy inequality, in that 
|n- 1 X?(k',z)tMX^(k, z)\ = (\/2){y/k'Jk z + yjk z /k' z ), which is greater than 1 
unless k' — k z . For closed channels the solutions are 



X((k;z) 



[m/ihKz)} 1 / 2 exp[-a(C)(wr/4 + K z z)\ 
a(()(l/2)(h 3 K z /m) 1 / 2 exp[-a(C)( i 7r/4 + k z z)] 



the inner products take the z-independent forms 

^- 1 ^'(k;z)tAf^(k;z) = 8? P 8 B < + 8? B 8 F <. 



(51) 



(52) 



In general, the properties that distinguish between between the four types 
of state of motion of a particle, that is open- versus closed-channel type, and 
F versus B type, depend on the local behavior of the state vector in wavenum- 
ber space (kt, k x , k y ). The corresponding position (t,x,y) space forms of these 
properties are nonlocal. As mentioned in Section 4, these properties are likely to 
complicate an attempt to make a physical interpretation, in the context of the 
present formalism, of measurements at a given z of local properties in position 
t, x, or y. This is in contrast to standard quantum mechanics with t as the 
evolution coordinate, where there is only one type of state in x, y, z: f -type 
and open channel. 

We note that the intermediate free-particle case 2mkt/h = k 2 
rise to a "dipole ghost" state, in that the reduced Hamiltonian on the rhs of 
Eq. 148(1 cannot be diagonalized by a similarity transformation. The construct, 
which Heisenberg named (see jSZ], references given therein, and jSH], P- 14), 
derives from the definition of the minimal polynomial of a finite-dimensional, 
square, complex matrix — see MacLane and Birkhoff 39 , Ch. IX. 6: Let I be a 
pseudo-Hermitean operator such that there exists a (real or nonreal) eigenvalue 
A of L, an integer n > 2, and a state X\ n so that the state (L — \) n ~ 1 X\ n 
is not the zero state and is an eigenstate in that [L — X) n X\ n — 0, then X\ n 



k 2 gives 
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will be called an "n-pole ghost" state of L associated with the eigenvalue A. 
We presume that, for any given L, and for each of its eigenvalues A, there is 
a bounded number — possibly zero — of types of ghosts associated with it. The 
ghost states associated with fixed L and A, and of different pole-orders n and m, 
are linearly independent of each other and of associated eigenstates; the direct 
sum of all the eigenstates and of all the corresponding linearly independent 
ghost states is a complete set of states in the overall space. 

Continuing with the zero-potential-energy, intermediate-case solutions, we 
note that the symbols F and B are not useful. We take the solutions X a (h(k 2 + 
k 2 )/(2m), k x , k y ; z) where a = 1, 2, as follows: 



X 1 {h{k 2 x + k 2 y )l{2m),k x ,k y -z) = 
X 2 (h(k 2 x + kl)/(2m),k x ,k y ;z) = 



[mp/K] 1 / 2 ' 


[m/(ph)]^ 2 2iz 

-i[fi 3 /Ml 1/2 



(53a) 
(53b) 



where p is an arbitrary positive number of dimension length introduced to make 
the components dimensionally consistent with Eqs. I|49l) and 151|) . The inner 
products are also z-independent: 

h-'X^mkl + k 2 y )/(2m) 1 k x ,k y ; z)^MX a {h(kl + k 2 y )/{2m) 1 k x ,k y ; z) 

fo, ifa' = a, (54) 
1+1, ifa'^a. 



Physically, the states in the intermediate case propagate parallel to any plane 
z=constant, that is, neither forward nor backward along z. Note that the so- 
lution Eq. I|53a|l is, and that of Eq. I|53b|) is not, an eigenstate with eigenvalue 



fact, the X 2 is a 



zero of the reduced Hamiltonian on the rhs of Eq. 
dipole ghost state for any choice of k x , k y and z. 

We next compute the inner product at each z of two free-particle wave 
functions ^(i, x, y, z) and x, y, z), when they have expansion amplitudes 
/^(k) and g^(k), respectively. It is convenient to divide k-space into domains 
for open and closed channels: 



d A k 



d A k = 



closed 



d A k = 



dk x dk y / dkt, 

Jh{kl+kl)/2m 
Mkl+k 2 y )/2m 

dk x dk y / dkt, 

J — oo 

d 3 k+ III d 3 k. 

open J J J closed 



(55a) 

(55b) 
(55c) 
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The inner product of W and $ is z-independent, and takes the form 



F 

(*(*);$(*)) = III (ft^O/'WVM 

JJJopcn <=B ^ 

d 3 k[f F (kyg B (k) + f B (ky g F (k)}. 

closed 

Note that the subspace generated by "i 7 "' open-channel states has a positive 
definite norm, while the space of U B" open-channel states has a negative defi- 
nite norm; each of these subspaces therefore comprises a Hilbert space. In the 
scattering phenomena analysed in Section 3 we shall discover that the open- 
channel sub-matrix of the 5*-matrix is unitary, and preserves the inner product 
of two vectors belonging to a direct sum of these Hilbert spaces referring to dif- 
ferent z-planes, assembled so that the sign of the inner product and metric are 
reversed in the second subspace component — hence there is a positive definite 
metric overall. The "F" and "B" Hilbert spaces on any z=constant plane are 
of limited utility, as most linear operators encountered in the space of states do 
not map such a Hilbert space into itself, but generate superpositions of F- and 
B-states, and of open- and closed-channel states. 

The question of normalizing the space-evolving wave functions can now be 
addressed: If the potential V(t, x,y, z) ^ in, and only in, the interior of the 
box £>2, a solution *&(t, x,y,z) of Eq. l]2.{[l can be expressed in an expansion of 
the type Eq. (|46|) . belonging to potential- free regions, in the neighborhood of 
both z = z\ and z = Z2, but with different expansion amplitudes at each end of 
the interval. We first define the basis functions 

E c (k;t,x,y,z) = (27r)~ 3/2 cxp(-ik t t + ik x x + ik y y)X c -{k; z). (57) 

We assume here and unless otherwise stated that there is no closed-channel 
input, and adopt the following conventions: 

*(t,x,y,zi) = III d 3 fc/£(k)S F (k;i,x,y,z 1 ) 

J J J open 

d 3 kC t (k)E B (k;t,x,y, Zl ), (58a) 
*(/..<■•//. .2 ) = /// d 3 kf i B n (k)E B (k;t,x,y,z 2 ) 

/open 

d 3 kf F At (k)E F (k;t,x 1 y,z 2 ), (58b) 

Since the flow of particle presence is conserved, the norms of ^(t, x, y, z\) and 
\l/(t, x, y, z 2 ) are equal: 

(!>*(«,) = (!>*(*!)• (59) 
Note that at z = z% , the forward and backward propagating parts of the wave 
function correspond to input and output, respectively, with the opposite asso- 
ciation at z — Z2- Note also that for a time-dependent potential energy there 
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will be scattering from open-channel input into both open- and closed-channel 
output, which circumstance is accounted for in Eq. 158|) . We now define a nor- 
malized, Type II wave function as one for which the input amplitude function 
is normalized to one, that is, 

1 = /// dMi/i^ + i/ifwi 2 ] 

JJJopen (g0) 

d 3 fe0/<fut(k)| 2 + |/ o B ut(k)| 2 ], 

/ open 

where the second equation follows from Eqs. H56[) . (|58|) . and (|59|l . Note that there 
is no contribution to the output normalization from closed-channel amplitudes. 

We conclude the section by developing a formula for the expectation value 
of the operator H2 in terms of the wave-number space expansion amplitudes 
/^(k). Let ^(t,x,y,z) be as in Eq. gSJ). Then we have 

{ih)^(z) = [[[ dtdxdy III d 3 k V x, y, z) f M 

J J Jr3 J J Jrs (61) 

x / c (k)(27r)~ 3/2 i exp[-ik t t + ik x x + ik y y]X c (k; z). 

Replacing t exp[— ik t t] by i(d/dk t ) exp[— iktt] and integrating by parts on k t , we 
find that 

F 

{thh(z) = ^ I / / dtdxdy d 3 k V *(t,x,y,z)tM 

JJJB 3 J J JR 3 £_g ^Q2~) 

1 d r 

x (2tt)- 3 / 2 exp[-iktt + ik x x + ik v y}-— / c (k)X^(k;; 

Analogous to spatial position operators in momentum space, the operator for t 
transforms into —id/dk tl the sign difference being a result of the negative sign 
in the exponent in Eq. H42|) . This result agrees with that in Ref. 001, Ch. 8, 
Eq. (286); see also [21], §8. Eq. 162fl can be interpreted as yielding a quantum- 
mechanical value for the average arrival time, or crossing time, of the particle 
in the state ^ at the given z=constant plane. 

The evaluation of Eq. I|62|l is facilitated by the following formulas: for £ = F 
and B and for Cf = B and F, respectively, and for open channels, 

rtC(k;z)tM 7W k;z) = (63a) 

ft- 1 X^(k;z)tAfi^(k;z) = <r(C)^ exp(- CT (C)2ifc^), (63b) 
while for closed channels 

n- 1 ^(k;z)tMi^(k;z) = - ff ( C )^ exp(-a(C)2^z), (64a) 

ft- 1 ^(k;z)tMi^i(k;z) = «r(0^. (64b) 

2 CrCf fiKz 
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If we carry out the differentiations in the integrand of Eq. I|62[l . we find that 



d 6 k 



open 



im 
2hk 2 
im 



exp(+2ifc z z)/ B (k)*/ F (k) 
exp(-2ik z z)f F (k)*f B (k) 



2hk 2 

+ S f (l/ F (k)| 2 + I/ B (k)| 



d A k 



closed 



r(k) 



a Of 1 



i dkt 



(k)+/ F (k 



B 



exp(-2 Kz z)|/ F (k)| 



i dkt 
m 

' 2Hk 2 



(k) 

cxp(+2 Kz z)|/ B (k)| 2 



^fB ikrf F {k) + !^l f F {kTf B {k) 

Tin? hn 7 



(65) 



Due to the denominators involving k 2 or k 2 in the above, convergence of the in- 
tegrals requires that the f' (k) approach zero sufficiently rapidly as k approaches 
the boundary between open and closed channels. 



3 Scattering; dwell and delay times 

In this section, we shall presume the presence of a generic potential energy 
distribution V(t,x,y, z), such that its support is contained within the box B2. 
The potential energy gives rise to scattering of the (we presume, purely open- 
channel) input signals, such that reflected and transmitted waves across the 
spectrum of k, including both open and closed channels, will comprise the out- 
put signal from the box. We now assume that our prescription for specifying 
the input yields necessary and sufficient information such that a solution to 
the Schrodinger equation within the box exists, satisfies the input boundary 
conditions, and is unique. Accordingly, the output is determined by the input, 
and this association must be linear in view of the linearity of the Schrodinger 
equation. The linear operator specifying this association consists of reflection 
and transmission coefficients, which can be assembled into an 5*-matrix, which 
in turn — as we shall verify — has a submatrix, referring to purely open-channel 
output as well as input, that is unitary. 

We presume that the Schrodinger equation has been solved for all open- 
channel inputs, and express the output linearly in terms of the input as follows: 
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/o S ut(k) 



/<fut(k) = 



d?k' 

open 

"^(k;k')/£(k') 
T BB (k;k')/^(k')' 

d 3 k' 

open 

T FF (k;k')/ F (k') 
fl FB (k;k')/£(k')' 



(66) 



(67) 



The functions T FF , R FB , R BF , and T BB arc reflection and transmission coeffi- 
cients, where the input-to-output superscripts are to be read from right to left. 
In Eqs. (|66|l and (|67|l . consistent with Eqs. I|58a|l and (|58b|l . the reflection and 
transmission coefficients are defined for the output parameter fc t having either 
an open- or a closed-channel value. 
For later convenience, we define 



/oP en (k / _ k „ } = 5 ope„ (fc / _ ^gfo _ _ ^ 



(68) 



where S opcn (k' t — fc") is defined only for both k' f and k' t ' corresponding to open 
channels. 

Let us now substitute Eqs. Ip3|) and l|B7|) into Eq. (JBDJ. We obtain a quadratic 
expression in the input amplitudes ff n on both sides of the resulting equa- 
tion. Since these amplitude functions are arbitrary, the coefficients of the four 
quadratic terms must be equal. We infer that, for both k' t and k' t ' being of 
open-channel type, 



d A k T Fi, (k;k')*T FF (k;k") 

open 



R a * (k;k')*iT^ (k; k") 
7 opcn (k' -k"), 



(69) 



d 6 k 



R FB {k; k')*i? FS (k; k") 

open 



T (k; k')*T BB (k; k") 
= 7 opon (k' -k"), 



(70) 



open 



R FB (k; k')*T FF (k; k" 



+ T SB (k; k')*i?^ (k; k") 



(71) 
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cPk 



R Bt (k; k')*T SB (k; k") 

open 



T^(k; k')*i?^(k; k") 



(72) 



0. 



Let us now make up an S'-matrix and its transpose conjugate from the 
reflection and transmission matrices. In the following, the unprimed index k t 
ranges over all real values, while k' t and k" range over open-channel values only: 



S(k;k") 
5t(k';k) 



T FF (k;k") R FB (k;k") 
R BF (k;k") T BB {k;k") 



(73) 
(74) 



~T FF (k;k')* R BF (k;k')*~ 
R FB (k;k')* T BB (k;k')*_ ' 

It is convenient to define two submatrices of S 1 , the open-channel part S and 
the closed-channel part S c as follows: 



S Q (k; k') = S(k; k'), for all k t > {h/2m){k 2 x + k 2 y ), 
S c (k; k') = S(k; k'), for all k t < (h/2m)(k 2 x + k 2 y ). 

Then S is unitary on the left as a result of Eqs. (|rJ9"|> - (|?2"|l : 

(SlS )(k';k")= fff d 3 k[Sl(k';k)S (k;k")] 

J J J open 

= h ®/ opcn (k' -k"). 
One expects that S is also unitary on the right, 

(S Sl){k!- k") = I 2 <g> 7°P cn (k' - k"). 



(75a) 
(75b) 



(76) 



(77) 



We now reduce the formulas for the expectation values of the operator il-i at 
z = z\,z%, using Eqs. (JHEJ), (E2J), JHEJ), and (frlTjl . and then establish a relatively 
simple form for the difference of the two values. We assume that the input am- 
plitudes ff n (k) and the output state values of the ^-matrix elements go to zero 
at the open/closed-channel threshold so that the following integrals converge. 
We have first 



open 



1 f)f F 

d»*/£(kr T -^(k) 

x/4(k)*^'(k;k';^)£(k0. 



d A k 



open 



(78) 



open 
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The matrix of coefficients is as follows: 
A FF (k;k'; Zl ) 



hk z 



r pon (k-k') 



+ ^pexp(-2 l fez 1 )i? BF (k;k') 
-i?^(k';k)*^exp(+^z 1 ) 
d 3 k" R BF (k";k)* 



R BF (k";k') 



open 

mz\ 1 d 

d 3 k" R BF (k";k)* 



closed 



X 2h^ eM2K " Zl)R (k " ;k ' ); 
A FB (k;k'; Zl ) 



im 



exp(-2ik z z 1 )T BB (k; k') 



/77 d 3 fc"i? BF (k";k)* 

J J J open 

TYLZ\ 1 9 



closed 



2?i< 2 



T BB (k";k') 
d 3 k" R BF (k"-k)* 
cxp(2<'z 1 )T siJ (k";k , ); 



^4 (k;k';^) 



= -T BS (k';k)^exp(2zfc^ 1 ) 
+ fff d 3 k"T BB {k"-k)* 

J J J open 

R BF {k"-k') 



open 

mzi 1 d 



+ fff d 3 k"T BB (k";k)* 

J J J closed 



m 



exp(2<si)ir* (k";k') 



(79) 



(80) 



(81) 
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A BB (k;k'; Zl ) 



-III 

J J J op 

mz\ 1 d 
III ^ 

J J J closed 



d 3 k"T BB (k";k)* 



T BB (k";k') 
"T BS (k";k)* 



2hn" 2 



exp^'zOT^k";^) 



At z — z 2 we have 



{ih)i$(z 2 ) 



fff d*kf B (V*l^ ; 

J J J op 



/ open 

d 3 k 



i dk t 

/open J J ./open ^ 

x/4(k)*C«'(k;k';^ 2 )/£(k'). 
In the above, the coefficient matrices are 
C FF (k;k';z 2 ) 

= Iff d 3 k"T FF (k";k)* 

J J J open 



mz2 1 d 



hk'l i dk'l 



closed 



T FF (k";k') 
d 3 k"T FF '(k";k)* 



TO 



2hn" 2 



cxp(-2K , z 'z 2 )T^(k";k'); 



C FB (k;k';* 2 ) 



= T FF (k';kr^eM-^Kz 2 ) 
fff d 3 k"T FF (k";k)* 

J J J open 

i? FS (k";k) 



open 

mz2 1 d 
hkj + ~iWl 



d 3 k"T FF '(k";k) 



closed 



TO 



2ft K " 2 



exp(-2<z 2 )i? FB (k";k') 
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C SF (k;k';* 2 ) 



2ftfc? 



exp(2ifez 2 )T FF (k; k') 



jff d 3 k" i? FS (k";k)* 

J J J open 



mz2 1 9 



T FF (k";k') 



/7Y d 3 /fc"i? FB (k";k)* 

J J J closed 

exp(-2<'z 2 )T FF (k";k'); 



2ft K " 2 



C BF (k;k';^ 2 ) 



mz 2 
hk z 
im 



7 opon (k - k') 
cxp(2ik z z 2 )R FB (k; k') 



+ i? FS (k';k)*^ 7 exp(-2 J ^ 2 ) 



/77 d 3 fc".R FB (k";k)* 

J J J open 

mz 2 1 9 



z dk'l 



closed 



i? FF (k";k') 
d 3 fc"i? FF (k";k')* 



■cxp(-2<^ 2 )i? FB (k";k'). 



2H K f 

The difference of the two expectation values takes the form 
(il 2 h (z2) -{il 2 ) Mzi ) = 1 1 1 d 3 k[[[ d 3 k' 

J J J open J J J open 

x^.^ n (k)^«'(k;k';z i; z 2 )/£(k'). 

C,C 

We break the D-matrices into constituents: 

D = Di + MD 2 S + S\D\M + & D 3 S + S\D A S C , 

where 

1 d mz\ 1 d mz 2 
hk' z ' i dk' t hk' z 



D\ = 7°P en (k-k')®diag(-^ 
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/ ^77 777 

D 2 = diagl ^2 exp(2ik z z 2 ), -^-^ exp(-2ik zZl ) ), MM) 
„ ,. fid mz 2 1 d mz\ , 



D 4 = diagf -^-^exp(-2<2 2 ),-^72 exp(2<2;i)J. (93) 

In Eqs. (EH) and the double primes indicate the dummy variables of inte- 
gration implicit in the final two summands on the rhs of Eq. 1)89(1 . The D\ and 
D3 terms were obtained by Smith , 021 for the case that the evolution coor- 
dinate is radial and the interaction is time-independent, but otherwise general. 
Note that Smith's 5-matrix is the transpose of the 5-matrix defined above. 

It is plausible that Eq. (|88() . given that the total input is normalized as in 
Eq. l|t)U|) , provides a complete expression for the mean dwell time of the particle 
in the box, inasmuch as it is also an expression for the space-time integral over 
the box of the divergence of the flow vector density of time. We remark that 
if the potential energy is time-independent, then the S^-matrix takes the form 
S (k";k') = S opcn (k' t ' - k' t )S (k' t ' ,k'^k^,k' x ,k' y ); one can now show that, due 
to the unitarity of S , the terms involving —idff n /dk' t (k.') cancel out in the 
overall expression for the dwell time. This cancellation does not, as we shall 
see, occur for the individual delay times for transmission or reflection from a 
zone of time-independent interaction. 

The average delay times that are measured in beam experiments for either 
transmission or reflection are not so fundamentally defined. We simplify the 
problem as follows: First, we assume that only one kind of input, that is F or but 
not and B, is present. Second, we assume that the closed-channel contributions 
will be negligible in the measuring apparatus. Third, we neglect interference 
between the incoming signal and the outgoing signal in the case of reflection 
(hence, the contributions linear in the S'-matrix are discarded) . We now specify 
what remains after these simplifications. 

In the first instance, let f£ (k) = 0, and let be normalized as in Eq. 
Pljl. The net out going reflected and transmitted currents of particle presence 
are called 1Z™?_ f (zi) and T F ^ F (z 2 ), respectively, and take the values 



rc3SL*(*i) = /// d 3 k / / / <Pk> / / / d 3 fc"[/£(k)* 

J J ./open J J ./open J J J open 

x (-l)R BF (k"; k)*i? BF (k"; k')/£(k')] , (94a) 
TrtM = 1(1 *kfff d*k> fff d 3 fc"[/ F (k)* 



x T FF (k"; k)*T FF (k"; k')/ F (k')] . (94b) 
According to Eq. jfiS} , we have 

T^ F (z 2 ) - K° B u i F (zi) = 1- (95) 
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The mean currents of time at entry and upon reflection at z\ , and upon trans- 
mission at z 2 , will be called, respectively, t f (zi), T B n ^ F (z\), and Tp$_ F {z%), and 
can be inferred from Eqs. (|79|) (twice) and (|84fl . subject to the three simplifica- 
tions spelled out in the previous paragraph, as follows: 



d 3 kf F (k) 



1 B^F 



open 

d 3 k 

open 

1 d mz\ 
idk? 4 



mz\ 



i_a 

i dkt ' hk 

F 



hk" 



^out / \ 
T F ^ F {Z 2 ) = 



(Fk 



open 



1 d mz2 

~iWl 4 



hk" 



open 

R BF (k";k') 



open 

FFf,ii,,l 



T**(k";k') 



/ta(k), (96a) 
d 3 k"R BF {k";k)* 

(96b) 
d 3 k"T FF (k";k)* 

(96c) 



open 



In the second instance, let /^(k) = 0, and let f B (k) be normalized as in Eq. 
<|6UI) . The net outgoing reflected and transmitted currents of particle presence 
are called Ti-° F ^_ B {z2) and T§^ B (zi), respectively, and take the values 



cPk 



d 3 k' 



d 3 k"[f B (ky 



i? FB (k";k)*i? FS (k";k')/^(k')], 
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According to Eq. (j7U|l , we have 

K7t B {z2) - nt B { Z1 ) 
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The net currents of time at z 2 upon entry and after reflection, and at z\ after 
transmission, will be called, respectively, t^(z 2 ), t f ^ b (z2), and t b u ^_ b (zi), and 
can be inferred from Eqs. (|87|) (twice) and (|82p . subject to the three simplifica- 
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tions given previously, as follows: 
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We now undertake to use the derived results to obtain estimates for the 
average delay time for the four processes of transmission and reflection. Due 
to the absence of space- and time-reversal symmetry of the potential energy, 
there will be no special relationships between the two transmission times or 
between the two reflection times. Let the transmission delay times be called 
r trans j- Z2 ^_ z ^ anc J Z2 ^ w hile the reflection delay times are called 

Tb^_ f (zi <— z\) and Tp^_ B (z2 <— 22)- We proceed from the following principle 
for computing delay times (currents are taken with their algebraic signs intact): 



delay time = 



(output current of time across exit plane) 
(output particle current across exit plane) 
(input current of time across entry plane) 



(100) 



(input particle current across entry plane) ' 
where the exit plane is the same, or the opposite, as the entry plane on reflection, 
or on transmission, respectively. (Similar formulas appear in |43| . p. 110, Eqs. 
(26) and (27), in Q2], Eqs. (28) and (29), in Eqs. (16) and (17), in 
Eq. (61), and in gS], Eq. (1).) We therefore have that 
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4 Discussion 

We have shown that the time can be construed as an observable in the context 
of the Schrodinger equation. We shall now undertake the discussion of some 
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aspects of the physical interpretation of the formalism, and of the relation of 
the present work to certain other similar investigations. 

For a state ^{t, x, y, z) made up of a finite number of closed-channel ampli- 
tudes that are all exponentially decreasing with z increasing, the net current of 
particle presence crossing any z=constant surface is zero — see Eq. ()56f) . Zero 
net current means that the particle will, in a large number of experiments, cross 
negatively as often as it crosses positively. The exponential decrease of the wave 
function is consistent with a picture of a "virtual" particle being created in the 
interaction region, propagating a small distance positively, and then falling back 
into the region where it was created. For F-type closed-channel states we say 
that the state propagates in the positive z direction, but represents an equal 
flux of particle motion in positive and negative z-dircctions. For open-channel 
-F-(-B-)type states, propagation and global particle motion are jointly positive 
(negative). 

It is possible, with the given physical system, to make a measurement of local 
(in (t, x, y)) net transit flux of particles across a proper subset of the z=constant 
plane. A local measurement of either the positive or the negative direction in 
z of a particle's transit involves noncommuting projection operators, one in po- 
sition space and one in wavenumber space. That is, we would then ask two 
incompatible yes/no questions: (1) Did the particle cross the plane in a given 
proper (t,x,y) subregion? (2) Did the particle cross positively/negatively in 
the z-direction? The outcome depends on the details of an often-repeated mea- 
surement on the system with the same input in each trial. To be sure, taking a 
sufficiently large subset of the z=constant plane for asking question (1) will, to 
a good approximation, be close to taking all of it, so that question (2) can be 
answered with negligible inconsistency. Also, the question "What is the differ- 
ence between the numbers of particles crossing positively and particles crossing 
negatively across a small subregion of the z-plane?" involves no inconsistency, 
and yields a result predictable from the local wave function alone; it is the sep- 
arate local densities of positive and of negative crossings that depend on the 
measurement scheme. 

A substantial effort has been dedicated to the establishment of a time-energy 
uncertainty principle — see the discussion and references in 07J . An uncertainty 
principle appears to be associated with a positive definite metric, a requirement 
that we have dropped. It is not obviously impossible to formulate some kind of 
analogous principle involving the time and its conjugate momentum pt within 
the present formalism in special circumstances, but the author has not been 
successful in finding such circumstances. 

The Schrodinger equation can create spatially and temporally localized ed- 
dies of probability current: even though a wave is made up entirely of a packet 
of F-type open-channel states, this current density can be negative in a neigh- 
borhood (see the examples by Kijowski and Mielnik cited in the following para- 
graph) and therefore will be greater than one in a complementary set of the 
given z=constant plane. Hence, even for a packet of only free-particle -F-type 
states this normalized particle count can fall outside the interval [0, 1] on proper 
subsets of a z=constant plane, and is not a probability, although a transition 
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process, as in the conventional interpretation of quantum mechanics, has ir- 
reducible randomness and is unpredictable in detail in contrast to a classical 
process. A further complication results from the circumstance that this normal- 
ized particle count is not, when closed channels are present, even globally (i.e., 
across an entire z=constant plane) an algebraic sum of two separate currents 
due to forward- and backward-propagating particles, as there can be nonzero 
global interference between the F- and B-type closed-channel contributions to 
the total current. 

Kijowski |24) . |48| undertook to establish a time-energy uncertainty princi- 
ple by analyzing the evolution of a Schrodinger wave function in a space-like 
direction, and in this respect there is overlap between Kijowski's work and the 
present undertaking. Kijowski's first "unsuccessful attempt" f[%4*|. §3) begins 
in a similar manner to that proposed above, but his inner product law does not 
involve an integral over time; since, as noted following Eq. H5()(l . the wave func- 
tions of Eq. (|49|l do not satisfy the Cauchy inequality, the interference terms in 
a local inner product can make the current density negative for the superposi- 
tion of two forward-propagating open channel states, as shown in an example 
in §3. Mielnik §5, Lemma) noted that a Schrodinger wave packet 

that at t = has its source entirely to the left of z = 0, say, could eventually 
give rise to probability currents normal to the z = plane that need not be 
everywhere positive. Similarly, the integrand for the particle current for the 
norm of a superposition of F-type open-channel states in Eq. I|28[l need not be 
everywhere nonnegative. These local negative currents all result from interfer- 
ence terms that yield zero net contribution in the present formalism due to the 
integral over t, x and y in the inner product. 

Kijowski's formalism is substantially different from the present one — the 
norm of an F-type state is given in [23] , Eq. (9) — but in which the average time 
of crossing a spatial wall for F-type states nevertheless reduces to the same form 
(EUi §10) m terms of the probability current as Eq. JSJJ. There are discussions 
of Kijowski's work in gH], §1.5.1, and |SD|, §10.2. 

Mielnik [20] critiques both Kijowski's [23] and Piron's [2S| attempts to es- 
tablish formalisms for spacewise propagation of a wave function, concludes that 
they do not offer a solution to the problem of defining time as an observable, 
and makes no additional proposals along these lines. Although the initial ideas 
of the two latter papers resemble that of the present paper, the respective imple- 
mentations differ considerably, so we shall not attempt further review of them 
here. 

Another question concerns the generalization of the effect of a measurement 
on a wave function that propagates in both directions away from the surface on 
which the measurement is performed. Suppose in fact that, in a problem of type 
II, two adjacent boxes occupy the space-like intervals [^1,^2] and [^2,^3], and 
that a measurement is made (over t, x, y) at z = z 2 , which measurement partly 
"collapses" the wave function there. The input at z-i to both boxes can change 
as a result of the acquired information, leading in turn to a change in the overall 
output at Z\ and Z3, and, due to reflections, a change in the wave function at z 2 
at which the measurement is made. There is therefore a problem of consistency, 



27 



in that the measurement at Z2 changes the outgoing wave function on both 
sides of Z2, and therefore, after a reflection, changes the ingoing wave function 
on both sides of z<z , that is, changes the frequency of results of the measurement 
at 22, and so on. This problem seems analogous to the "grandfather paradox" 
(see jS^, Ch. 4) of the influence of a physical system with itself between two 
different i=constant surfaces, when reflection as well as transmission of signals 
along the evolution coordinate occurs. No successful attempt at analysis of this 
class of measurement problems is known to the author. 

We reemphasize that it is the Schrodinger equation that is taken as funda- 
mental in the present argument, and that a theory of measurement, a probability 
interpretation, and an uncertainty principle, are all presumed to be derivative 
ideas that may require alterations from their conventional forms in order to 
bring them into concord with the body of formalism presented here. We have 
assumed that not just ip*ip, but also the spatial components of the probability 
four-current of Sec. 2, are measurable quantities. Moreover, we assume that the 
flux density of a physical quantity represented by an operator to is given by the 
four-current of Eqs. In classical terms, the total fluxes amount to normal- 

ized counts of particles crossing a given oriented three-dimensional surface in 
space-time, weighted by the time of crossing or some other such quantity, and 
also weighted positively or negatively according as the particle crosses positively 
or negatively when it transits the given surface. Note that in experimental tri- 
als the particle is presumed always to transit positively across a segment of a 
£=constant surface — Galilean geometry singles out inconstant planes from all 
other planes in space-time and permits special treatment for these cases. 

To recapitulate in other words, we assert that the evolution of a Schrodinger 
wave function in a spatial direction does not generally admit of description 
in terms of probability amplitudes. The claim is rather that the wave func- 
tion in type II problems permits only the computation of certain expectation 
values, that is, average results of many repeated experiments with the same 
input signal, but such that there exists no underlying distribution of nonncga- 
tive quantities, analogous to ip*t[>, that accounts for the results. We advocate 
the non-introduction of the term "negative probabilities" as there seems to re- 
sult a decrease in physical clarity thereby; instead, a kind of random behavior 
more general than that which can be characterized by probabilities is entailed. 
Khrennikov, in j52], Ch. III. 2 and references given therein, describes what he 
calls "signed 'probabilistic' measures", which could serve as a classical analog 
of the local norm of Eq. Ij26(l : see also [23, §34. We infer that the conven- 
tional, probability interpretation of nonrelativistic quantum mechanics should 
be subordinated to an interpretation involving observable stochastic currents 
of particle presence, and, more comprehensively, observable stochastic currents 
of other physical quantities as temporal position, spatial position, energy, mo- 
mentum, and so on. This "particle current" interpretation of the formalism can 
describe systems of both types I and II; the usual probability interpretation 
then applies in problems of type I and other special cases. 

Although what appears to be a mathematically consistent formalism has 
been constructed herein, and a preliminary physical interpretation advanced, 
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many questions along these lines need to be addressed, and consistency with 
experimental tests established, before the proposal can with confidence be re- 
garded as a physical theory. 

The above limitations notwithstanding, the formalism proposed herein has 
obtained results that agree to an extent with some special results previously 
derived, and has secured results that would be difficult to obtain by other pub- 
lished methods of analysis: for example, a generic expression Eqs. (|9~3l) 
for the average dwell time for a particle reflecting from or passing through a 
time-dependent barrier. 
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